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We construct a non-Grassmann spinning-particle model which, by analogy with quantum me- 
chanics, does not admit the notion of a trajectory within the position space. The pseudo-classical 
character of the model allows us to avoid the inconsistencies arising in the quantum-mechanical 
interpretation of a one-particle sector of the Dirac equation. 
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Non-abelian gauge groups play a crucial role in field 
theory as well as in the Standard model. In this work 
we observe that they imply new possibilities when used 
in the construction of finite-dimensional theories as well. 
We suggest and discuss the notion of pseudo-classical me- 
chanics (pCM), the term by which we refer to models 
with a number of observable configuration-space vari- 
ables less than the number of physical degrees of free- 
dom. In other words, we consider classical-mechanical 
models which, by analogy with quantum mechanics, do 
not admit the notion of a trajectory within the position 
(i.e. configuration) space. So we expect that pCM turns 
out to be useful in describing quantum phenomena by 
(semi) classical methods [1-11]. Classical mechanics with 
such a strange property can be constructed on the basis 
of a singular Lagrangian with a multi-parametric group 
of local symmetries. As examples of pCM we present 
the models which are invariant under transformation of 
the non-abelian gauge group with two and three local 
parameters. Symmetries imply functional ambiguity in 
solutions to equations of motion: besides the integration 
constants a, the solution depends on the arbitrary func- 
tions eair), X — f{T,Ci,ea{T)). According to the general 
theory of singular systems [12-14], variables with ambigu- 
ous dynamics do not represent observable quantities. So, 
when we are dealing with the locally-invariant theory, our 
first task is to find candidates for observables, which are 
variables with unambiguous dynamics. Equivalently, we 
can look for the gauge-invariant variables. 

We start in Section 2 with a couple of toy models and 
show that, generally, it is impossible to construct the ob- 
servables within the position variables only. It is worth 
noting that on the phase space there always is the well- 
defined notion of a trajectory [13]. In Section 3 we con- 
sider a more realistic case, presenting the non-Grassmann 
model of the Dirac electron. In Section 4 we show how 
the pseudo-classical character of the model allows us to 
solve the problems arising [2, 3, 18] when we try to ap- 
ply the methods of relativistic quantum mechanics to a 
one-particle sector of the Dirac equation. 



One of the local symmetries which will be presented in 
our models is reparametrization invariance. So, we first 
outline the reparametrization invariant formulation of a 
relativistic particle. 

The motion of a particle in special-relativity theory 
can be described starting from the three-dimensional ac- 
tion —mcj dtJc^ — (^)^- The problem here is that 
the Lorentz transformations, x'^ = A'^j^x", act on the 
physical dynamical variables x*(t) in a higher nonlinear 
way. To improve this, we pass from three-dimensional 
to four-dimensional formulation. Introducing the para- 
metric representation x^{t) = {ct{T) , x'^ (t)) for the tra- 
jectory x'(i), the particle can be described by the La- 
grangian action 



4e^ ' 2 ' 



S= I dri^an'' - ^m'c'). (1) 

The corresponding Hamiltonian action reads 



Sh = drpf.i'' + Pee - -e{p^ + rr?c^) - X^Pe, (2) 

where Ae(T) stands for the Lagrangian multiplier of the 
primary constraint tTe = 0. Variation of the functional 
implies the Hamiltonian equations 



A. 



Pe 



0, 



ep'' 



Pm = 0, (3) 



as well as the constraints iVe = 0, p"^ + rn^c^ — 0. We 
note that the variable Ae(r) cannot be determined with 
the constraints, nor with the dynamical equations. As 
a consequence (see the first of Eqs. ([3])), the variable 
e turns out to be an arbitrary function as well. Since 
e(r) enters into the equation for x^ , its general solution 
contains, besides the arbitrary integration constants, the 
arbitrary function e(r). Hence the only unambiguous 
ones among the initial variables are p^ and tTe, see Eqs. 
([3|). x'^ has one-parameter ambiguity due to e. 

The ambiguity reflects the freedom in the choice of 
parametrization for the particle trajectory {a(T) is an 
infinitesimal function) 
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T ^ T = T — a, 
2:^(r) ^.T'^(r') = x''(t), 
e(r) ^e'(T') = (l + d)e(T). 



then 



5e 



ax^ 



(ae)'. 



(4) 



Action ([T]) turns out to be invariant under the 
reparametrizations. 

By construction, the expression for the physical trajec- 
tory x'(i) is obtained resolving the equation x° — x^(t^ 
with respect to r, r = t{^^^i then x'(t) = x*(t(x°)). 
Using the expression 



d! /(r) 



dx" x^(t) ' 



(5) 



for the derivative of a function given in parametric form, 
we obtain 



dx' 
~di 



y/jf + m^c 



2„2 



(6) 



Eq. ([6]) coincides with that of the three-dimensional for- 
mulation. As should be the case, the physical coordinate 
x*(i) has unambiguous evolution. 

Toy model which admits the position-space tra- 
jectories. Consider the Lagrangian action 



5 = 



2ei 



(7) 



This is written on the configuration space x^^, e\ and 62; 
the Minkowski metric is r\^^ = (— , -t-, -|-, -t-). The action 
is invariant under the reparametrizations, 5x^ = ai^, 
6ei = (aei)', Se2 = (Q!e2)', as well as under the following 
transformations with the parameter /3(r) 



Sx'^^Px'', 5ei^2f3ei, 



(562 - fi- 



(8) 



The transformations form a non-abelian group, [(5q, Sp] — 
5n, P = —a/3. Each local symmetry removes two degrees 
of freedom [2^ , so the number of configuration-space ob- 
servables is equal to 2. 

The Hamiltonian of the theory ([7|) is 



H 



ei 2 



e2{px). 



(9) 



This implies the Hamiltonian equations (in what follows, 
we omit equations for the auxiliary variables e^, Pefc, as 
they are not necessary for discussion of the ambiguity of 
the variable x^ ) 



as well as the constraints 

p2 ^ 0, {px) = 0. 



-62?'^, 



(10) 



(11) 



The equation for x^ has two-parametric ambiguity due 
to ei and 62, while that for p^ has one-parametric ambi- 
guity. Inside the light-cone, we construct the variables 



\J —x^ 



r 



Their equations read 



X 



ei 



V- 



~2' 



V^ 



p 



x'^p^ . 



0. 



(12) 



(13) 



They can be compared with Eqs. ([3]). The constraints 
([Tl]) acquire the form p^ = 0, {px) = 0. Note that x is 
/? -invariant variable. So the ambiguity presented in Eq. 
(|13p is due to the reparametrization symmetry. In ac- 
cordance with this observation, we assume that the func- 
tions x^{t) and p'^''{t) represent the reparametrization- 
invariant variables x'^{t) andp'^(i) in the parametric form. 
Their equations of motion read 



dt 



pu 



-pO' 



dpi" 

~dr 



= 0. 



(14) 



Since they are unambiguous, the variables x^{t) and p^{t) 
are candidates for the observables. 



By construction, x^ obey the identity 5:^5:^ 



-1. So 



only three of them can be taken as coordinates of the 
configuration space. Adding the variable a = .^ ,^ to 



the set x', we obtain a coordinate system. As the two 
independent observables we can take the gauge-invariant 
variables x^it) and a;^(i). Hence, the present model ad- 
mits observable trajectories constructed within the posi- 
tion space; see Eq. ((T2|) . 

To conclude, we point out that Poisson brackets of 
the Lorentz-covariant observables ([T^ generate the non- 
commutative algebra 

{i^,p''}=iV^''(i), {x^,i''} = 0, {p^,p'^}=p[^i''l(15) 

Here and below we denote 



^^"^(a) = r/^"- 



(16) 



Toy model without position-space trajectories. 

Consider the following Lagrangian action written for the 
variables x'', lo^ , ei and 62 



S ■ 



dr 



1 



2(ei-e2) 



[{Dxf + 2e2{DxLj) + exC? 



-m c H — uj 
2 2 



(17) 



We have denoted Dx^^ 



62^'' 



The non-abelian 



gauge group is composed by reparametrizations as well 
as by the following transformations with the parameter 

/3(r): 



(5x" = 



/? 



,{Dx^' + e2Lo^'), de^ ^ jS. (18) 



ei 

This implies that the number of physical degrees of free- 
dom on configuration (phase) space is equal to 6 (12). 

Denoting conjugate momenta of x, lo by p, n, the 
Hamiltonian of the theory ([TT]) reads 

i/ = Itt^ - ic.2 + iei(p2 + m'c^) + e2p^(c.^ - 7r^)(19) 

This implies the Hamiltonian equations 

x" = eip" + e2{uj'' - tt"), p" ^ 0, 



w'^ = vr^ - e2p'', t:'' =uj'' - e2p''. 



as well as the first-class constraints 



P 



2 2 

m c 



0, 



{P,^- 



0. 



(20) 



(21) 



The equation for x has two-parametric ambiguity due to 
ei and 62, while those for uj and tt have one-parametric 
ambiguity. 

Taking into account the first-class constraints, we 
could expect 6 observable dynamical variables on the 
configuration space. However, it is easy to see that 
any configuration-space quantity a^{x,uj) with one- 
parametric ambiguity is proportional to w''. Similarly 
to the previous model, this can be used to construct 
only three unambiguous dynamic variables. As the six- 
dimensional configuration space can not be spanned with 
the unambiguous variables, the model represents an ex- 
ample of pseudo-classical mechanics. 

On the phase space we can construct various variables 
with one-parametric ambiguity due to 62 



(px) 



-P' 



X 



62 (W^ 



(22) 



{i^^Q'') 



rt^n)^ 



pt^Oj" 



(pn) p2 



{J^^c:>"} = r, 






{J^'•',n°'}^'n"^''n''l 



(pit) 



(27) 



The set of 12 independent observables of the phase- 
space can be selected as follows. We parameterize the ini- 
tial space by the coordinates x°, x*, p'', w'', 5' = e^^'^Jjk 
and 7 — \/{uJ + 7r)2. The dynamics of the theory is re- 
stricted on the surface pp' + m^ = 0, (pui) = which is 
invariant under the action of the gauge group^. The sur- 
face can be parameterized by x'^, i*, p*, w', S*' and 7. 
The corresponding dynamic variables i'(i), p^{t), w*(t) 
and S^{t) have unambiguous dynamics. Hence we can 
take them as the independent observables. 



III. NON-GRASSMANN MECHANICAL 
MODEL OF THE DIRAC ELECTRON 



P^ 



/ = 0, 



(w + 7r,p)^ 



2p^ 



{u^' - tt'^). 



U! 



(23) 



^e^iu^-TT^), (24) 



^i^JTW 



TT = .J . (25) 



y(Z7T^ 



J^"' ^LU^'^^'' -uj^wt', jt"" ^ e2P^^'iuj - tt^I (26) 



The constraints 1^1^ acquire the form p^ + m?c^ — 0, 
{pCj) = 0. The new variables are invariants of the /3- 
transformation^. So the ambiguity presented in Eqs. 
(|22p - (j26p is due to the reparametrization symmetry. Sim- 
ilarly to the case of the relativistic particle, we assume 
that the functions x'^(r), p'^(r), w''(t), ^/^(r) and J^"'{T) 
represent the physical variables x*(i), p^{t),... in the 
parametric form. According to Eq. ([5]), the dynamics 
of the physical variables is unambiguous. 

The Poisson-bracket algebra of the Lorentz-covariant 
observables is highly noncommutative; the nonvanishing 
brackets are 



{^^x'^} 



p[pi^] 



{w^w'^} 



p[M^H 






As a more realistic example of pCM, we discuss the 
spinning-particle model suggested in a recent work [15]. 
The configuration space of the model consist of the dy- 
namical variables Q"(t) — {x^^ , uj^ , ui^) as well as the 
auxiliary variables e;, Z = 1, 2, 3, 4. x'^ are coordinates of 
the Minkowski space with the metric 77'*" = (— , +, +, +). 
The spin-space uj^ = (cj'*, w^) is equipped with 5*0(2, 3) - 
metric rj^^ — (—,+,+,+,—). Consider the Poincare- 
invariant Lagrangian 



1 



-GapQ'^Q 



64 A 

2 ^ 



6/ 



-ai. 



(28) 



We have denoted ai = m c ,a2 = mch, and 03, 04 are 
real numbers. In what follows, we discuss the free theory. 
Interaction with an external electromagnetic field will be 
discussed at the end of Section 4. The kinetic term looks 
like that of a free particle moving on the curved nine- 
dimensional space with the metric 



Gap = I -e2W^G^i/ eiG^i, 

62 



esA 



-7^v 



63 A 



■iOy 



6?W" 

B 

63 A 



(29) 



We have denoted G^^ = ^V\p.v - Sai;^^], 
e\{iJ'f and >1 = B -K e\{uj^'f. 
We introduce the abbreviation 

Dx^ = x^ - — (w^w^ - w^w^), 

63 



B 



6163 



(30) 



^ Due to the identities i'^p^ = 0, Tr'^Tfp 
not all of them are independent. 



1 and e'^'^Jij Jofc 



^ We use the phase-space form of reparainetrizations, &ui^ 
—ae2{uj'^ — tt''), 5p>^ = 0; see [15] for details. 



then the Lagrangian ((28)) can be written as follows, 



63 



2e3 



"24 4, 



ei 



^-UJ-UJA - -^aiiSl) 



The Lagrangian is invariant under a three-parametric 
group of local symmetries. One of them is the 
reparametrization symmetry. Besides, there are two 
more symmetries with the local parameters /3(t), 7(t) 

S^x'^^Pp'^, 5^ei=/3; (32) 



S^u! — je^n , (5-y7r = —7640; , 
^763 = (763)', ^764 = iiei)'. (33) 

Here p^ = -3^ , tta = -Mx ■ In the Hamiltonian formula- 
tion, the Lagrangian (|3T)) leads to the following Hamilto- 
nian [15] 



H 

£3 
2 



ei 



62, 



(p2 + m^c^) + ^{Pf.J^'' + mch)+ 

64, 

2 



(tT^TTa + 03) + —{UJ'^UJA + 04) + AeaTTea, (34) 



UJ 



e^TT^ +e2{puj), 



e2{pTr) 



03 5 

04 



(39) 



The only unambiguous variable is p'^. The configuration 
space variables x^, ui have two-parametric ambiguity. 

Let us compute the total number of physical degrees 
of freedom. Omitting the auxiliary variables and the cor- 
responding constraints, we have 18 phase-space variables 
x^^, pp, w^, TTA subject to the constraints (|36)) - (|38)) . Tak- 
ing into account that each second-class constraint rules 
out one variable, whereas each first-class constraint rules 
out two variables, the number of physical degrees of free- 
dom on the phase space is 18— (2-1-2 x 3) — 10. Hence we 
could expect five observables on the configuration space. 
However, using the configuration-space variables only, it 
is impossible to construct five unambiguous quantities. 
Thus, once again we have an example of pseudo-classical 
mechanics. 

Let us discuss the physical sector of the phase space. 
A brief inspection of the equations of motion allows us 
to construct the Lorentz-covariant variables with one- 
parameter ambiguity. They are the five-dimensional 
angular-momentum tensor*^ J'^^ — 2a;['^7r-^l; this obeys 



where 



J^" ^2{uj^t:''~lj''t:^). 



(35) 



If we omit the spin-space coordinates, cu — tta — 0, the 
Hamiltonian reduces to that of the spinless particle, see 



The Hamiltonian implies the constraints 



iU OJa + 04=0, 



TT^UA = 0; 



(36) 



j5t-^-e2 J^^Pu , j^'" = -€2 f" [^''1 ; 



as well as the position variable 



2p'^ 



X — ep^ . 



(40) 



(41) 



We have denoted e = ei + |^- So the reparametrization- 
invariant variable x^ (t) has a deterministic evolution: 



p^ + vr?(? = 0, TT^TTA -I- as = 0; 



P^J' 



mch, = 0. 



(37) 



(38) 



The first one states that configuration space of spin 
is anti-de Sitter space. The constraints ([36]) form the 
second-class pair while those of Eqs. (|37| and (|38|) are 
the first-class constraints. The constraint (|38l) . being im- 
posed on the state vector, leads to the Dirac equation^, 
(7''Pp + "^c)^ = (see [15] for details). 

The Hamiltonian equations of the theory read 



1 



eip^ 



762 



J^^, 



p^ 



w'^ = eaTT^ + 630; V, T^^ 



= 0; 

04 



dx^ 
~di 



pj 



(42) 



As the classical four-dimensional spin vector, we take the 
Pauli-Lubanski vector which has no precession in the free 
theory 



gt. 



1 



fiuaP 



Pu^a 



/3' 



S^" = 



(43) 



In the rest frame p^ — (mc, 0,0, 0), it reduces to the 
three-dimensional rotation generator, 5° = 0, 5' = 
■^mce^^^Sjk, as is expected in the non-relativistic limit. 

We point out that the second term in Eq. (|^T|) has the 
structure typical for non-commutative extensions of the 
usual mechanics, see (2311. 



The present model implies both the Dirac equation and the mass- 
shell condition p^ + m^c^ = 0. The model without the mass- 
shell condition has been discussed in [16, 17]. This shows the 
same undesirable properties as those of the Dirac equation in 
the classical limit [2, 3]. 



* Note that the constraints H36|l and H37|l fix the value of the 
Casimir operators of the 50(2, 3) group. Besides, they guar- 
antee that J^f^ is the time-like vector (J^*^)^ = —4(03(0;^)^ -|- 
a^ln^)^) < 0, for positive values of 03, 04; see [15] for details. 



IV. CONCLUSION. PSEUDO-CLASSICAL 

MECHANICS AND THE CLASSICAL LIMIT OF 

THE DIRAC EQUATION 

Although a true understanding of spin is achieved 
in the framework of quantum electrodynamics, a lot of 
effort has been spent in attempts to construct a me- 
chanical model of a spinning electron, see [1-11, 15- 
17] and references therein. The Dirac spinor ^ can 
be used to construct the four-dimensional current vec- 
tor, ^7^^, which preserves for solutions to the Dirac 
equation, 9^ (^7^^) ~ 0. Hence its null-component, 
^'^^ > 0, admits the probabilistic interpretation, and 
we expect that a one-particle sector of the Dirac equa- 
tion could be described in the framework of relativistic 
quantum mechanics (RQM). 

However, it is well known that adopting the RQM in- 
terpretation, we arrive at a rather strange and controver- 
sial picture [2, 3, 18]. To recall this to mind, we use the 
Dirac matrices a* and /3, to represent the Dirac equation 
in the form of the Schrodinger one 



ihdt'^ = i7*, H = ca'pj + TOc^/3. 



(44) 



Then H may be interpreted as the Hamiltonian. If we 
pass from the Schrodinger to the Heisenberg picture, the 
time derivative of an operator a is ihd = [a, H], For the 
basic operators of the Dirac theory we obtain 



ihdi — 2{cpi — Hai 



K = 0. (45) 



Below we enumerate the inconsistences arising in the 
RQM interpretation of these equations and show how 
our pCM allow us to avoid them. 

The wrong balance of the number of degrees of freedom. 
Assuming x as the position operator, the first equation 
in (1451) implies that the operator cd^ represents the ve- 
locity of the particle. Then the physical meaning of the 
operator p^ became rather obscure in both the semiclas- 
sical and the RQM framework. Using the quantum-field- 
theory arguments, Foldy and Wouthuysen [6] argued that 
the basic operator x which appears in the Dirac equation 
does not correspond to the observable quantity. Then 
they constructed the position operator X* with reason- 
able properties. 

Our model supports the Foldy- Wouthuysen suggestion. 
Indeed, we observe that the variable a;^ is not a gauge- 
invariant quantity in our model, so it is not an observable. 
The observable variable which can be associated with the 
particle position is x^{t). According to Eq. (|42t . p* de- 
termines its velocity. x^{t) written in Eq. (j4ip repre- 
sents the Lorentz-covariant analog of the operator X* in 
the classical theory. We also point out that the Foldy- 
Wouthuysen transition a; — >■ X corresponds in pCM to 
the transition from the gauge-non-invariant to the gauge- 
invariant variable. 

Zitterbewegung. The equations (|45p can be solved, 
with the result for a;'(i) being [2, 3] x' = a* + bp^t + 



c*exp(— 2i^i)_ The last term on the r.h.s. of this equa- 
tion states that the free electron experiences rapid oscil- 
lations with higher frequency ^ ~ ^^1^- ^^ is often as- 
sumed that Zitterbewegung represents the physically ob- 
servable motion of a real particle [18]. The analogous 
systems that are described by a Dirac-type equation and 
simulate Zitterbewegung are under intensive study in dif- 
ferent physical set-ups, including graphene, trapped ions, 
photonic lattices and ultracold atoms, see [i9| and refer- 
ences therein. 

Our model excludes the Zitterbewegung, since this rep- 
resents the dynamics of the unobservable variable x^. 
The observable variable x* moves along a straight line, 
see Eq. (|42l) . 

Velocity of an electron. Since the velocity operator ca* 
has eigenvalues ±c, we conclude that a measurement of 
a component of the velocity of a free electron is certain 
to lead to the result ±c. 

In our model, the conjugate momentum p^^ determines 
the velocity of the physical coordinate S^, see Eqs. ((4T|) . 
(1321) • Then the mass-shell condition (|57)) guarantees that 
the particle cannot exceed the speed of light. 

Bargmann- Michel- Telegdi (BMT) equations. In their 
seminal work [7], Bargmann, Michel and Telegdi sug- 
gested relativistic equations for the classical trajectories 
and spin precession in uniform fields. The equations 
practically exactly reproduced the spin dynamics of po- 
larized beams and agreed with the calculations based on 
the Dirac theory. While the BMT model apparently does 
not imply the Dirac equation, there is a certain relation- 
ship between the two schemes. Namely, the first term of 
the WKB solution to the Dirac equation can be used to 
construct the quantities which obey the BMT equations 

MM- 

We have seen above that performing the canonical 
quantization of our model in the initial variables we ar- 
rive at the Dirac equation. Now we show that physical 
variables of the model obey the BMT equations. 

We take the Hamiltonian of interacting theory in the 
form 



H 



ei 



{V' + —F^uJ'"" + m'c')+ 



62 



2c 



{Vf,J^^ + mchy 



— -(tT TTA -1-03) -I- — -(w WA + 04) + AeaTTe 



(46) 



where Vu 



'-.A^j^ is the mechanical momentum. 



^ — P/i 

This does not break the local symmetries presented in the 
model for the case of uniform electric and magnetic fields. 
However, the interaction deforms the physical sector: un- 
ambiguous variables of the free theory no longer remain 
unambiguous in the interacting theory. In particular, 
x^ , p^ and S^ now have two-parametric ambiguity. Up 
to the order 0{h ), the quantities with one-parametric 
ambiguity turn out to be 



■pA" 



2cV'- 



■{FJV), 



-eF^.P" 









S^ = — eF'^^S''. (47) 



We have denoted e = ei — fpje2. The correspond- 
ing reparametrization-invariant variables obey the BMT 
equations with g = 2 

(j P^ d e ■ 

_x' - r— —P' - —F'-'^'P 

dt " PO' dt " PO "" 



^S^ = -^P'"^'^- (48) 
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